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We present results for the mass of the flavor singlet meson calculated on two-flavor full QCD conflgurations 
generated by the CP-PACS full QCD project. We also investigate topological charge fluctuations and their 
dependence on the sea quark mass. 



1. INTRODUCTION 

Recently considerable progress has been made 
in the simulation of full QCD ||l[. In particu- 
lar, sea quark effects have been found to lead to 
a closer agreement of the light meson spectrum 
with experiment |^,^. 

Missing from the calculated spectrum, how- 
ever, has been the flavor singlet meson ?/. Due 
to the difficulty of the determination of the dis- 
connected contribution, only preliminary lattice 
results have been available ■ 

In the first half of this article, we present new 
results on this problem. Since these are obtained 
with two flavors of dynamical quark, we call the 
flavor singlet meson as rj and reserve the name 77' 
for the case of Nf~3. 

The rj' meson is expected to obtain a large mass 
through the connection to instantons. This leads 
us to an investigation of topology in full QCD, 
presented in the latter half of this article. 

Calculations have been performed on configu- 
rations of the CP-PACS full QCD project |||. 
These have been generated using an RG- improved 
gauge action and a tadpole-improved SW clover 
quark action at four different lattice spacings and 
four values of the sea quark mass corresponding 



Table 1 

Parameters of lattices used in this calculation. 



N'^xNt 


P 


TVs a [fm] 


a ^ [GeV] mps/my 


12^x24 


1.8 


2.58(3) 


0.92(1) 


0.81-0.55 


16^x32 


1.95 


2.45(3) 


1.29(2) 


0.80-0.59 


243 X 48 


2.1 


2.60(5) 


1.82(3) 


0.80-0.58 


243 X 48 


2.2 


2.06(6) 


2.30(7) 


0.80-0.63 



to mps/mv~ 0.8-0.6. An overview of the simula- 
tion parameters is given in Table |^. More details 
can be found in [|,|. 

2. FLAVOR SINGLET MESON 

The mass difference Am between the fiavor sin- 
glet (77) and non-singlet (tt) meson can be ex- 
tracted from the ratio 



m = 



disc 



1 - Scxp(-ATOi), (1) 



*talk presented by R. Burkhalter at Lattice '99 



where the right hand side indicates the expected 
behavior at large time separation t. 

In this work, the connected propagator was cal- 
culated with the standard method. For the dis- 
connected propagator we used two methods. In 
the first instance, it was calculated using a vol- 
ume source without gauge fixing (the Kuramashi 
method) g. This measurement was made after 
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Figure 1. Ratio for /3 — 1.95 and mTJirip — 0.59. 
The solid line is from a fit with Eq. n. 

every trajectory in the course of configuration 
generation for all runs listed in Table |l|. As for 
the second method, we employed a U(l) volume 
noise source with 10 random noise ensembles for 
each color and spin combination. This was per- 
formed only at /3 = 1.95 on stored configurations 
separated by 10 HMC trajectories. 

Figure |l] compares the ratio R{t) for the two 
methods. We observe that they are consistent 
with each other but that the error is smaller for 
the first method. This might be due to the fact 
that there are 10 times more measurements with 
it, although binning is made over 50 HMC tra- 
jectories for both cases to take into account auto- 
correlations. In the following we only use data 
obtained with the first method. 

In Fig. [I] we also see that the error of R(t) in- 
creases exponentially, which make the determi- 
nation of Am via Eq. |l| impossible at large time 
separations. The data, however, shows the ex- 
pected behavior already beginning from small t 
and a fit with Eq. ^ is possible from t„iin — 2. 
Increasing tmin leads to stable results, as can be 
seen in Fig. ||. 

The chiral extrapolation of mf^ linear in the 
quark mass is shown in Fig. ^. Contrary to the 
pion mass, the flavor singlet meson remains mas- 
sive in the chiral limit. 

Figure ^ shows the rj meson mass at all mea- 
sured lattice spacings after the chiral extrapola- 
tion. The scale is set using the p meson mass. 
A linear extrapolation to the continuum limit 
gives rrijj — 863(86) MeV. This value lies between 
the experimental 77(547) and 77' (958) masses. We 
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Figure 2. Chiral extrapolation of 77 meson mass. 
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Figure 3. Continuum extrapolation of 77 mass. 

emphasize that a proper comparison with ex- 
periment requires the introduction of a third 
(strange) quark and a mixing analysis. 

3. TOPOLOGY 

Studies of topology on the lattice have encoun- 
tered several difficulties. In addition to the am- 
biguity of defining a lattice topological charge, 
it was found that topological modes have a very 
long auto-correlation time in the case of full QCD 
with the Kogut-Sussind quark action 0. 

We employ the field theoretic definition of the 
topological charge together with cooling. For 
the charge we use a tree-level improved definition 
which includes a lx2-plaquette, hence the O(a^) 
terms are removed for instanton configurations. 
For the cooling we compare two choices of im- 
proved actions, both including a 1 x 2-plaquette 
term: 1) a tree- level Symanzik improved (LW) 
action and 2) the RG improved Iwasaki action. 

Using different actions for cooling can lead to 



3 



N,=0 



' cool with LW action 
cool with RG action 



0.00 0.05 0.10 0.15 0.20 

Figure 4. Topological susceptibility for Nf = 0. 

different values of the topological charge. This 
ambiguity is only expected to vanish when the 
lattice is fine enough. We have tested this explic- 
itly by simulating the pure SU(3) gauge theory at 
three lattice spacings in the range a ~ 0.2-0.1 fm 
and at a constant size of 1.5 fm. As Fig. ^ shows, 
the topological susceptibilities Xt — {Q'^)/V for 
the two cooling actions converge to a common 
value towards the continuum limit. Using ^/a ^ 
440 MeV we obtain xt = (178(9) MeV)^, in agree- 



ment with previous studies 

In full QCD we have so far measured the topo- 
logical charge at /3 = 1.95. Figure ^ shows 
the time history for two quark masses. Auto- 
correlation times are visibly small even for the 
smallest quark mass. For the Wilson quark ac- 
tion rather short auto-correlation times have been 
reported in Ref. The fact that we find even 
shorter auto-correlations might be explained by 
the coarseness of our lattice. 

Based on the anomalous flavor-singlet axial 
vector current Ward identity, one expects the 
topological susceptibility to vanish in the chiral 
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Figure 6. Topological susceptibility at [3 ■ 

limit. Indeed, Fig. || shows the width to be shrink- 
ing with the quark mass. The decrease, however, 
is not sufflcient; as we find in Fig. |^, the dimen- 
sionless ratio Xt /c^ , with a calculated for each sea 
quark mass, does not vary much with the quark 
mass, and takes a value similar to that for pure 
gauge theory. To understand the origin of this 
behavior, more investigations at different lattice 
spacings will be needed. 
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Figure 5. Charge history at (3 = 1.95. 



